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ABSTRACT 

-V-  The  establishment  of  a  quasi-steady  flow  in  a  generic  scramjet  combustor  i3  studied  for 
the  case  of  a  time  varying  inflow  to  the  combustor.  Such  transient  flow  is  characteristic  of 
the  reflected-shock  tunnel  and  expansion- tube  test  facilities.  Several  numerical  simulations 
of  hypervelocity  flow  through  a  straight-duct  combustor  with  either  a  side-wall-step  fuel 
injector  or  a  centrally-located  strut  injector  are  presented.  Comparisons  are  made  between 
impulsively  started  but  otherwise  constant  flow  conditions  (typical  of  the  expansion- tube 
or  tailored  operation  of  the  reflected-shock  tunnel)  and  the  relaxing  flow  produced  by  the 
“uudertailored”  operation  of  the  reflected-shock  tunnel.  Generally  the  inviscid  flow  features, 
such  as  the  shock  nit  tern  and  pressure  distribution,  were  imaffected  by  the  time  varying  inlet 
conditions  and  approached  steady  state  in  approximately  the  times  indicated  by  experimental 
correlations.  However,  viscous  features,  such  as  heat  transfer  and  skin  friction,  were  altered 
by  the  relaxing  inlet  flow'  conditions,  A-'  y  i,  .  .-i  A  -  Cor.  i  ;  Corr  ° 
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1.  INTRODUCTION 


The  National  Aero-Snae.e  Plane  i'NASP'I  nroiect  lias  been  the  focus  of  the  recent  revival 

i  V,  ■  y  i  J 

of  hypersonic  aerodynamics  and  propulsion  research.  Propulsion  studies  for  the  high  Mach 
number  range  ( Mfught  >  4)  have  concentrated  on  the  air-breathing  supersonic-combustion 
ramjet  (scramjet)  using  hydrogen  as  the  primary  fuel.  Most  of  the  early  research  studies 
of  the  mixing  and  combustion  processes  in  a  scramjet  relied  on  conventional  ground  test 
facilities  in  which  the  test  gas  is  heated  to  flight  enthalpy  by  combustion  of  hydrogen  or  with 
an  electric  arc.  These  tests  have  generally  been  limited  to  Mjught  <  8  because  of  the  limited 
total  temperatures  available  with  these  facilities.  Generation  of  test  flows  with  enthalpies 
relevant  to  the  flight  regime,  10  <  Mfught  <  25,  requires  a  different  approach. 

One  type  of  facility  capable  of  this  high  enthalpy  range  is  the  pulse-type  facility  in  which 
the  test  gas  is  rapidly  heated  by  the  passage  of  a  strong  shock  wave.  Two  examples  of 
pulse  facilities  are  the  reflected-shock  tunnel  and  the  expansion  tube.  A  review  of  these 
facilities  and  their  use  in  experimental  hypervelocity  aerodynamics  has  been  presented  by 
Stalker  [1],  When  using  pulse  facilities  to  simulate  high  enthalpy  flows,  a  number  of  trade¬ 
offs  between  test  time  and  test  gas  conditions  must  be  made.  One  of  the  difficulties  in  using 
pulse  facilities  to  test  supersonic  combustors  is  that  the  test  times  are  relatively  short  -  on 
the  order  of  one-half  millisecond.  These  short  test  times  are  a  concern  because,  even  in 
steady  flow,  many  of  the  important  flow  processes  in  the  combustor  will  take  significant  time 
to  become  established  (i.e.,  approach  steady  state).  With  a  transient  flow  the  establishment 
may  be  further  delayed. 

Much  effort  has  been  directed  toward  the  optimization  of  the  scramjet  combustor  and, 
in  particular,  to  the  efficiency  of  fuel-air  mixing  and  reaction  so  as  to  obtain  minimum 
length  and  wall  heat  load  in  the  engine.  Because  of  the  high  velocity,  high  temperature  and 
low  density  of  the  combustor  flow,  the  mixing  of  the  fu  d  and  air  is  relatively  slow  at  the 
higher  flight  speeds.  Consequently,  mixing  augmentation  and  enhancement  by  the  controlled 
generation  of  turbulence  and/or  large-scale  vortices  in  the  mixing  and  flame- holding  zones  are 
being  examined  in  both  numerical  studies  [2]  and  experimental  studies  [3].  These  techniques 


generally  involve  some  type  of  recirculating  flow  in  the  wake  of  a  body,  step  or  transverse 
jet  and  hence  require  some  finite  time  to  reach  steady  state  (or  a  quasi-periodic  state). 

The  focus  of  this  study  is  to  examine  several  time  accurate  simulations  of  steady  and 
transient  flow  through  a  generic  scramjet  in  an  attempt  to  assess  the  following  issues: 

9  The  time  required  for  the  flow  in  the  wake  region  of  the  fuel  injector  (without  injection) 


4o  become  established  relative  to  the  typical  test  times  available  in  pulse  facilities 
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•  The  impact  of  test  flow  relaxation  (as  encountered  in  the  undertailored  operation  of  a 
reflected-shock  tunnel)  on  the  approach  to  steady  state  combustor  flow. 

•  The  validity  of  using  the  “hypersonic  equivalence  principle”  (see,  e.g.,  [4])  to  transform 
the  relaxing  flow  to  a  quasi-steady  flow. 


The  flow  through  the  supersonic  combustor  was  computed  at  the  high  enthalpy,  short 
duration,  and  sometimes  transient  conditions  typical  of  a  shock  tunnel  or  expansion  tube. 
Calculations  of  the  combustor  flow  were  m..de  in  a  time -accurate  manner  using  the  Navier- 
Stokes  code  SPARK  [5],  [6].  This  code  has  been  previously  applied  to  the  simulation  of 
scramjet  flows  with  fuel  injection  and  a  thrust  ramp  in  a  shock  tunnel  environment  [7].  In 
the  present  study  the  code  is  used  to  compute  the  basic  combustor  flow  without  fuel  injection. 


2.  OPERATION  OF  PULSE  FACILITIES 

A  brief  description  of  the  operation  of  two  particular  pulse  facilities  will  be  given  in  this 
section.  The  first  is  a  free-piston  driven  reflected-shock  tunnel,  T4,  located  at  the  University 
of  Queensland  [8],  The  second  is  an  expansion  tube,  HYPULSE,  located  at  General  Applied 
Science  Laboratories,  New  York  [9]. 


2.1.  Reflected-Shock  Tunnel 


The  principal  features  of  a  shock  tunnel  and  its  operation  are  shown  in  Figure  1.  The 
driver  tube,  which  initially  contains  low  pressure  helium,  and  the  shock  tube  which  contains 
the  test  gas  are  separated  by  the  primary  diaphragm.  Attached  to  the  downstream  end  of 
the  shocktube  is  a  nozzle  whose  throat  is  significantly  smaller  than  the  diameter  of  the  shock 
tube.  The  test  gas  is  contained  by  a  thin  mylar  diaphragm  which  separates  the  shock  tube 
from  the  evacuated  test  region. 

The  first  stage  of  operation  of  the  T4  facility  is  the  launch  of  the  piston  and  its  acceleration 
down  the  compression  tube.  The  driving  force  is  supplied  by  the  compressed  air  from  the 
reservoir.  The  primary  diaphragm,  which  is  typically  composed  of  two  sheets  of  2mm  mild 
steel,  subsequently  bursts  at  a  pressure  56.6ik/Po.  At  this  point,  the  helium  has  been 
compressed  to  1  / 60th  of  its  initial  volume  and  is  contained  in  approximately  0.5m  of  the 
compression  tube. 

After  diaphragm  rupture,  helium  driver  gas  expands  into  the  shock-tube  and  compresses 
the  low  pressure  test  gas  before  it.  The  incident  shock  travels  the  length  of  the  shock 
tube,  reflects  from  the  closed  end  of  the  tube  and  brings  the  test  gas  to  rest  in  the  nozzle 
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supply  region.  Operation  in  this  manner  is  called  “tailored”  [10]  and  is  shown  in  the  (x- 
t)  diagram  (Figure  1(b))  by  the  contact  surface  coming  to  rest  when  intercepted  by  the 
reflected  shock.  The  compressed  test  gas  is  contained  in  the  nozzle  supply  region  in  a  length 
of  approximately  0.25  —  0.5m.  Ideally,  the  nozzle  supply  conditions  are  maintained  as  the 
reflected  shock  continues  upstream  through  the  driver  gas. 

Upon  shock  reflection,  the  light  secondary  diaphragm  bursts  and  some  of  the  compressed 
test  gas  expands  into  the  nozzle.  Smith  [11]  has  provided  a  quasi  one-dimensional  model  for 
the  nozzle  starting  processes. 

For  combustion  studies,  the  useful  test  time  is  terminated  by  the  contamination  of  the  test 
gas  by  the  driver  gas.  The  mechanism  for  this  contamination  is  the  bifurcation  of  the  reflected 
shock  into  two  weaker  oblique  shocks  near  the  wall  of  the  tube.  This  mechanism  has  been 
studied  in  [12]  and  has  been  shown  to  occur  for  undertailored  conditions  at  high  enthalpies 
[13].  For  undertailored  operation,  conditions  are  such  that,  when  the  reflected  shock  reaches 
the  helium-air  interface,  it  accelerates  into  the  helium  and  an  expansion  propagates  into  the 
nozzle  supply  gas.  With  no  other  influences,  the  supply  region  will  reach  a  new  equilibrium 
but  there  will  be  a  significant  and  unavoidable  drop  in  nozzle  supply  pressure  P,  shortly 
after  shock  reflection.  The  expansion  of  the  test  gas  slug  delays  the  arrival  of  helium  jetting 
down  the  walls  of  the  shock  tube  and  prolongs  the  useful  test  time.  The  approximate  time 
at  which  contamination  of  the  test  gas  is  expected  to  reach  the  nozzle  throat  decreases  with 
increasing  stagnation  enthalpy,  H,,  (see,  e.g.,  [1]). 

The  effect  of  finite  driver  size  upon  the  relaxation  of  P,  is  also  noticeable.  Stalker  [14] 
has  suggested  a  mode  of  operation  where  the  motion  of  the  piston  (after  rupture)  maintains 
approximately  constant  conditions  in  the  driver  tube  while  nearly  50%  of  the  helium  flows 
into  the  shock  tube.  Following  this  initial  period,  of  approximately  1ms,  the  pressure  in  the 
driver  decays  rapidly  and  this  decay  is  transmitted  downstream  along  a  (u-f  o)  characteristic. 
The  effect  of  the  finite  volume  of  the  driver  tube  is  that  the  the  decay  in  PM  may  start  early 
in  the  test  period. 

The  net  result  is  that  the  history  of  the  nozzle  supply  pressure  and  a  representative  pitot 
pressure  at  the  exit  plane  of  the  Mach  5  nozzle  will  appear  as  shown  in  Figure  2.  The 
pressure  transducer  measuring  Pt  is  located  a  few  centimetres  upstream  of  the  reflecting  end 
of  the  shock  tube.  Its  trace  shows  the  passage  of  the  incident  shock  and  the  arrival  of  the 
reflected  shock.  Due  to  transducer  rise  time  and  location,  the  full  reflected  pressure  is  not 
recorded.  Past  the  peak  value,  there  is  a  continuous  decay  in  P,  due  to  the  combined  effects 
of  undertailored  operation  and  driver  pressure  decay.  For  helium  driver  gas  and  operation 
considered  here,  this  decay  is  typically  25-30%  during  a  nominal  0.5ms  test  time  [15].  Note 
the  delay  between  the  supply  pressure  trace  and  the  pitot  pressure  trace.  This  delay  will  be 


3 


used  in  the  normalization  procedure  discussed  in  Section  3.5. 


2.2.  Expansion  Tube 


An  expansion  tube  facility  consisting  of  a  driver  tube  (initially  containing  high  pressure 
driver  gas),  a  shock  tube  (containing  the  test  gas)  and  an  acceleration  tube  (containing  low 
pressure  acceleration  gas)  is  shown  schematically  in  Figure  3.  The  driver  tube  and  shock 
tube  are  separated  by  the  primary  diaphragm  and  the  acceleration  tube  is  separated  from 
the  shock  tube  by  a  secondary  (and  very  light)  diaphragm.  The  shock-  and  acceleration 
tubes  have  the  same  diameter.  Although  the  HYPULSE  facility  has  a  fixed  driver  tube 
rather  than  a  free-piston  driver,  the  gas-dynamic  processes  in  the  shock  tube  are  similar. 

The  operation  is  initiated  by  the  rupture  of  the  primary  diaphragm.  High  pressure  driver 
gas  flows  into  the  shock  tube  and  the  incident  shock  compresses  and  accelerates  the  test  gas. 
On  reaching  the  end  of  the  shock  tube,  the  secondary  diaphragm  bursts  (ideally  without 
causing  disturbance)  and  the  shock- compressed  test  gas  expands  into  the  the  low  pressure 
gas  in  the  acceleration  tube.  The  pressure  to  which  the  acceleration  gas  rises  is  below  the 
pressure  of  the  bulk  of  the  shock-compressed  test  gas  and  so  the  downstream  portion  of  the 
test  gas  undergoes  an  unsteady  expansion  to  the  test  flow  conditions. 

The  test  time  commences  after  the  test-gas/acceieration-gas  interface  arrives  and  usually 
finishes  with  the  arrival  of  the  downstream  end  of  the  unsteady  expansion  fan.  Figure  4 
shows  a  typical  history  of  the  static  pressure  in  the  test  region.  The  trace  shows  a  rapid  rise 
due  to  the  shock  through  the  acceleration  gas,  a  slower  rise  (possibly  due  to  diffuser  starting 
processes)  to  the  test  conditions  at  approximately  60 fis  and  then  roughly  constant  until  the 
arrival  of  the  unsteady  expansion  400/i.s  later. 


3. 


NUMERICAL  SIMULATION 


The  computations  of  flow  through  the  generic  supersonic  combustor  were  made  with  the 
Langley  Research  Center  SPARK-2D  code.  This  two-dimensional,  elliptic,  finite  difference 
code  was  developed  by  J.  P.  Drummond  [5],  [6]  to  integrate  the  conservation  equations  for 
mass,  momentum  and  energy  in  a  time  accurate  manner.  Since  the  primary  concern  is  the 
temporal  development  of  the  combustor  flow  features,  the  scope  of  the  modelling  was  reduced 
by  considering  only  laminar  flow  of  a  nonreacting  test  gas  and  avoiding  the  fuel  injection  and 
mixing  issues.  The  effects  of  transition,  turbulence,  mixing  and  chemistry  will  be  examined 
in  future  studies. 
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3.1.  Combustor  Geometry  and  Boundary  Conditions 

A  generic  scramjet  combustor,  typical  of  those  tested  in  the  reflected  shock  tunnel  T4 
[15],  is  shown  in  Figure  5.  This  model  consists  of  the  major  features  of  a  scramjet  combus¬ 
tor,  including  some  form  of  injector  (with  an  associated  wake/mixing  region)  and  a  set  ol 
confining  walls  (with  their  associated  boundary  layers). 

To  perform  the  two-dimensional  computations,  the  scramjet  combustor  duct  in  Figure  5 
was  modelled  as  shown  in  Figure  6.  Only  half  of  the  duct  was  considered  in  the  computational 
domain.  Two  classes  of  calculations  were  performed  which  differed  in  the  location  of  the 
inflow  boundary.  Cases  1,2,3  and  6  with  the  inflow  boundary  at  the  start  of  the  duct  walls 
(x  =  0),  have  a  flow  domain  as  shown  at  the  top  of  Figure  6.  These  are  called  the  “short” 
injector  cases.  Cases  4  and  5,  with  the  inflow'  boundary  at  the  injector  strut  leading  edge 
(x  =  —0.138m),  are  called  the  “long”  injector  cases  (Figure  6(b)).  Descriptions  of  the  six 
cases  considered  in  this  study  are  shown  in  Table  1. 

Boundary  conditions  on  the  wall  (at  y  —  0.02357m)  and  the  surfaces  of  the  injector  were 
set  to  zero  velocity,  zero  normal  pressure  gradient  and  a  fixed  temperature  T  =  300if  (i.e., 
cold  wall).  Cases  4  and  5  included  a  free  boundary  along  y  =  0.02357m  for  —0.138m  <  x  < 
0m.  Conditions  at  the  inflow  plane  were  supersonic.  The  outflow  boundary  at  x  =  0.5m 
was  specified  as  supersonic  with  zero  gradients  in  the  flow  direction.  In  Cases  1  and  2,  the 
x-axis  downstream  of  the  injector  (x  >  0.114)  was  set  as  a  cold  wall  so  as  to  provide  a 
situation  close  to  that  found  in  parallel  wall  injection  studies.  Gases  3,4,5  and  6  set  this  line 
to  be  a  symmetry  plane  with  zero  normal  velocity  and  zero  normal  gradients  for  pressme, 
temperature,  and  tangential  velocity. 

3.2.  Grids 

The  grids  used  in  the  computation  are  shown  in  Figure  7.  The  entire  domains,  including 
the  injectors,  were  included  in  the  grids  and  during  the  computation  all  grid  points  within 
the  injector  were  reset  to  initial  conditions.  Each  grid  consisted  of  a  number  of  separate 
zones  in  which  the  nodes  were  clustered  at  one  or  more  of  the  zonal  edges  so  that  large  flow 
gradients  near  the  boundaries  and  planes  defining  the  injector  could  be  resolved  without 
excessive  computational  expense. 

The  x  and  y  coordinates  for  the  grids  used  in  the  short  injector  cases  were  independently 
generated  using  one  of  Roberts’  [16]  stretching  transformations.  (See  also  [17],  Section  5-b.l). 
In  each  zone  there  were  N  +  1  nodes,  including  the  end  points,  located  at 

z  =  ZaT)  +  2{,(1  -  77),  (1) 
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where 


=  [(/?  +  2a)A  -0  +  2a] 

(2) 

V  ~  (2a -i  1)  [1  +  A]  ’ 

4  i\("-“)/(1-) 

(3) 

V  =  =  0...JV. 

(4) 

Details  of  the  zonal  boundaries  and  stretching  parameters  are  given  in  Table  2.  These 
parameters  were  adjusted  until  numerical  artifacts  were  eliminated  from  flow  monitors  such 
as  the  length  of  the  recirculation  region  behind  the  injector  (discussed  in  Section  4.3)  and 
the  position  of  the  recompression  shock  (discussed  in  Section  4.4).  Note  that,  for  a  =  0,  the 
nodes  are  clustered  close  to  the  za  end  of  the  zone  and,  for  a  =  0.5,  the  nodes  are  clustered 
close  to  both  boundaries. 

For  the  long  injector  cases,  the  y-coordinates  were  generated  in  a  similar  manner.  How¬ 
ever,  the  x-coordinates  were  generated  using  a  compression  function  due  to  Thomas  et  al. 
[18].  In  each  2one  there  were  N  +  1  nodes  (including  end  points)  specified  by  Equation  1 
and 

a  _  sinhifa-Aq)  ,  J  /eN 

’  =  '  + 1  ■ 


+  1  . 


where 


A  1,  r  l-f  (e^3- 1)770  1 

n  2  _1  +  (e~0  -  1)7?0 

and  Tj0  is  the  value  of  7/  at  which  maximum  compression  occurs.  Details  of  the  zonal  bound¬ 
aries  and  stretching  parameters  are  given  in  Table  3. 


3.3.  Initial  Conditions  and  Inflow  Conditions 

The  initial  conditions  within  the  combustor  were  set  to  T  —  3U0 K ,  P  —  3Q0.Fa  and  zero 
velocity  so  as  to  approximate  the  conditions  in  T4  experiments  in  which  the  test  section  of 
the  facility  is  evacuated  to  approximately  2  Torr. 

The  inflow  conditions  were  set  to  approximate  the  flow  in  the  T4  facility  fitted  with  a 
Mach  5  contoured  nozzle.  The  nozzle  supply  pressure  was  approximately  P,  =  52MPa  and 
total  enthalpy  was,  H,  =  SAM  J/ kg.  Details  of  the  experiment  are  recorded  in  [15]. 

The  test  flow  conditions,  both  in  the  shock- reflection  region  and  the  nozzle,  were  com¬ 
puted  as  quasi-steady  flows.  The  nozzle  supply  conditions  (behind  the  reflected  shock)  were 
estimated  with  the  FORTRAN  program  “ESTC”  [19]  which  incorporated  an  equilibrium 
model  for  air.  From  the  shock  reflection  conditions,  the  gas  expanded  adiabatically  (and  in 
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chemical  equilibrium)  to  the  measured  nozzle  supply  pressure  of  52 MPa.  Using  these  condi¬ 
tions,  the  flow  at  the  exit  plane  of  the  nozzle  was  estimated  with  the  quasi-one-dimensional 
program  “NENZF”  [20]  in  which  the  test  gas  consisted  of  a  mixture  of  the  species  N2,  N, 
02l  O,  NO ,  NO+  and  e~ .  The  gas  was  assumed  to  be  in  chemical  equilibrium  at  the  nozzle 
throat  but  a  finite  rate  chemistry  mode!  with  11  reactions  was  used  in  the  expansion  region 
of  the  nozzle. 

The  “steady”  conditions  at  the  end  of  the  nozzle  (and  inlet  to  the  combustor)  were  spec¬ 
ified  as  free-stream  velocity  3670ms-1,  a  static  pressure  of  21.5 kPa,  a  static  temperature  of 
1165jT,  and  oxygen  and  nitrogen  mass  fractions  of  0.2314  and  0.7686  respectively.  Although 
the  NENZF  program  determined  a  species  mixture  that  included  finite  amounts  of  0  and 
NO,  only  a  nonreacting  mixture  of  02  and  N2  was  used  in  the  SPARK-2D  computations. 

The  “transient”  inflow  was  constructed  from  the  nozzle  supply  pressure  trace  shown  in 
Figure  3,  assuming  a  quasi-steady  flow  and  using  the  ESTC  and  NENZF  codes  to  compute 
the  nozzle  exit  conditions  for  a  number  of  values  of  Pe.  The  time  evolution  of  P„  was  then 
approximated  by  the  curve  fit 

Pm  =  Pp«,fc(0.215  +  4.540f  -  22.148*3),  0  <  t  <  0.28, 

=  Ppeak(exp(-0.77H{t  -  0.28))),  0.28  <  t  <  2.26,  (7) 

where  P^k  =  62 MPa,  t  is  in  ms  and  P„  is  in  M Pa.  The  conditions  at  the  inlet  to  the 
scramjet  combustor  were  then  related  to  P,  with  the  following  equations 


Pxn[kPa)  =  0.435^-1.55, 

(8) 

Tin(K)  =  525  Pj-256, 

(9) 

f/m(ms-1)  =  350Tm°333. 

(10) 

The  histories  of  the  inlet  conditions  are  shown  in  Figure  8. 
3.4.  Numerical  Damping 


Some  difficulty  was  experienced  in  starting  the  solution  with  such  low  pressure  and  large 
incident-shock  Mach  number  (which  was  approximately  10  for  Case  1).  The  difficulty  ap¬ 
peared  to  be  caused  by  the  int  Taction  of  the  strong  incident  shock  structure  with  the  newly 
forming  recompression  shock.  To  get  the  calculation  started,  the  CFL  number  was  linearly 
varied  from  0.1  to  0.8  over  the  first  18000  time  steps.  For  Cases  2,3,4  and  5,  the  coefficient 
for  the  artificial  viscosity  terms  was  set  to  a  constant  value  of  -1.0  while,  for  Cases  1  and  6,  it, 
was  set  to  -1.2.  To  faithfully  capture  heat  transfer,  temperature  smearing  was  not  included. 
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3.5.  Normalization  Procedure  for  Transient  Flow 


The  liow  quantities,  q[x,  t),  computed  in  the  transient  Cases  2,  3  &  5  were  normalized  with 
upstream  reference  values  to  produce  equivalent  quasi-steady  values,  q( x,  t ).  The  particular 
transformation  used,  based  on  the  “hypersonic  equivalence  principle”  was 


q(x,i)  = 


q(x,t) 


(11) 


0)  t  Tjelay') 

where  the  time  delay,  rjelau,  was  given  as  the  nominal  transit  time  from  the  reference  position 


^dclay  — 


X  ~  Zo 
£/.„  ' 


(12) 


Essentially  this  means  that  individual  parcels  of  fluid  were  followed  downstream  through  the 
combustor  and  measured  quantities  were  referenced  to  the  parcels’  initial  states.  (See,  e.g., 
Section  4.8  of  [4].) 


4.  RESULTS  AND  DISCUSSION 


In  previous  studies  of  impulsively  started  flows  [21]  [22]  [23]  [24]  one  principal  result 
has  been  the  time  taken  for  the  flow  to  approach  steady  state.  This  establishment  time 
re  is  usually  combined  with  a  characteristic  velocity,  UCI  and  a  length,  Lc,  to  form  the 
dimensionless  parameter 

C  =  IjT  (13) 

which  represents  the  ratio  of  the  flow  establishment  time  to  the  time  for  flow  to  proceed 
through  the  domain.  For  the  flow  features  in  the  generic  scramjet,  such  as  boundary  layers 
and  recirculation  regions,  G  will  have  different  values.  For  the  impulsively  started  flow  over 
a  flat  plate,  where  the  boundary  layer  is  laminar  and  Lc  is  the  distance  from  the  start  of 
the  plate,  G  is  approximately  3  [21].  When  tin-  boundary  layer  is  turbulent  (7  ~  2.  For  the 
wake  of  a  sphere,  with  Lc  as  the  length  of  the  recirculation  region,  Holden  [22]  gives  G  ~  30, 
based  on  the  pressure  measurements  and  G  ~  70  based  on  heat  transfer  measurements. 

In  this  section,  six  numerical  simulations  of  the  transient  flow  in  a  generic  scramjet  com¬ 
bustor  will  be  examined  with  emphasis  on  the  approach  to  steady  state  (or  quasi-steady  state 
as  determined  by  the  hypersonic  equivalence  principle).  The  following  parameters  will  be 
used  to  determine  flow  establishment  times:  the  boundary  layer  thicknesses  at  the  end  of  the 
injector  strut,  the  length  of  the  recirculation  region,  the  location  of  the  recompression  shock 
and,  the  pressure,  heat  flux  and  shear  stress  on  the  y  =  0.02357m  wall.  The  computation 
of  the  flow  development  for  each  of  the  cases  was  continued  to  a  time  of  approximately  1ms 
which  is  typical  of  the  flow  time  in  the  current  generation  of  pulsed-flow  facilities. 
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4.1.  Overview  of  the  Solutions 

Figures  9  an'1  10  show  the  pressure  fields  for  the  short  injector  geometry  with  steady 
inflow  (Case  1)  and  transient  inflow  (Case  2)  respectively.  Contour  spacing  is  2 kPa,  These 
figures  highlight  the  similarity  between  the  steady  inflow  and  the  transient  inflow  solutions. 
The  pressure  distributions  along  the  y  =  0  and  y  =  0.02357m  walls  are  shown  on  the  right 
side  of  each  figure.  Note  that  the  stretching  of  the  y-coordinate  (by  a  factor  of  5)  changes 
the  appearance  of  the  shock  and  Mach  angles. 

At  the  start  of  the  duct,  there  are  a  pair  of  weak  shocks  generated  by  the  initially  rapid 
growth  of  the  boundary  layer  on  the  duct  walls.  Although  not  evident  in  the  pressure  field 
plots,  boundary  layers  develop  along  the  upper  wall  and  the  injector  surface.  The  boundary 
layer  on  the  injector  separates  at  the  step  (x  =■-  0.114m)  and  reattaches  to  the  y  —  0 
boundary  some  distance  downstream.  The  expansion  from  the  corner  of  the  injector,  and 
the  recompression  shock  formed  near  the  end  of  the  recirculation  zone  are  clearly  visible  for 
both  cases.  Subsequent  reflections  of  the  recompression  shock  from  the  duct  walls  can  be 
seen  for  t  >  0.2ms. 

The  initial  shock  structure  sweeping  through  the  duct  for  t  <  0.16ms  is  composed  of  an 
incident  shock  propagating  downstream  through  the  quiescent  gas,  a  contact  discontinuity 
separating  the  initial  duct  gas  and  the  incoming  test  gas,  and  an  upstream-facing  shock 
decelerating  the  incoming  test  gas.  Both  the  second -order  and  fourth-order  MacCormack 
schemes  experienced  some  difficulty  with  this  strong  shock  structure  and  numerical  oscilla¬ 
tions  are  evident  in  the  solutions  for  t  <  0.20ms.  However,  it.  is  expected  that  the  solutions 
are  reliable  at  later  times. 

The  incident-shock  position  was  identified  as  the  point  at.  which  the  wall  pressure  rose 
from  (the  initial  value  of)  300Pa  to  800Pa.  The  position  was  recorded  every  100 /is  and  the 
data  points  were  then  fitted  with  the  straight  line 

x  =  v(t  -  t0).  (14) 

The  values  of  v  and  t0  are  given  in  Table  4.  In  all  cases,  the  shock  speed  was  virtually 
constant  throughout  the  flow  domain.  However,  the  transient  inflow  cases  (2,3  and  5)  had 
a  lower  incident-shock  speed  than  the  steady  inflow  cases.  For  Cases  1  and  6,  the  position 
of  the  upstream-facing  shock  was  identified  as  the  point  where  the  wall  pressure  crossed  a 
35 kPa  threshold.  The  velocities  for  this  shock  are  also  given  in  Table  4. 

As  shown  by  the  pressure  contour  plots  and  the  wall  pressure  plots  in  Figures  9  and 
10,  the  overall  flow  features  of  both  steady  inflow  and  transient  inflow  cases  are  similar. 
Although  the  starting  shock  structure  took  longer  to  exit  the  duct  for  the  transient  inflow 
(Case  2),  the  major  steady-state  features  such  as  the  boundary-layer  interaction  shocks,  the 


9 


expansion  fan  and  the  rccomprcssion  shock  are  in  virtually  the  same  places.  The  principal 
difference  seems  to  be  the  average  level  of  the  duct  pressure  which  can  be  seen  relaxing  in 
the  wall  pressure  plots  of  Figure  10. 

The  effect  of  including  the  complete  injector  strut  in  the  computational  model  can  be 
seen  in  Figure  11,  which  shows  a  few  frames  from  Cases  3  and  5.  The  frames  are  placed 
back  to  back  and  the  times  have  been  chosen  to  show  approximately  the  same  stages  of 
development.  However,  the  instantaneous  inflow  conditions  are  not  exactly  the  same  in  each 
pair  of  frames.  The  major  differences  are  the  strong  shock  and  expansion  propagating  from 
the  front  of  the  long  injector.  Also,  the  shock  propagating  from  the  leading  edge  of  the 
y  =  0.02357m  wall  is  stronger  than  the  boundary-layer  interaction  shocks  seen  in  the  short 
injector  simulations.  After  reflecting  from  the  injector  strut  and  then  the  duct  wall,  this 
leading-edge  shock  interacts  with  the  expansion  fan  propagating  from  the  base  of  the  strut. 

Since  there  is  only  a  small  difference  between  the  long-injector  and  short-injector  simula¬ 
tions,  the  following  discussion  will  concentrate  on  the  details  of  the  short-injector  simulations. 

4.2.  Injector-Strut  Boundary  Layer 

For  the  short  injector  simulations,  the  boundary  layer  along  the  strut  surface  closely 
approximates  the  ideai  flat-piate  boundary  layer,  the  main  difference  being  the  boundary- 
layer  interaction  shock  striking  the  strut  surface  ati  ~  0.09m.  This  weak  shock  did  not  seem 
to  affect  the  flow  profiles  within  the  boundary  layer  at  the  end  of  the  strut  (x  =  0.114m) 
where  the  following  thicknesses  were  calculated: 

•  total  boundary  layer  thickness  8;  the  y- distance  from  the  strut  surface  to  0.99  of  the 
local  free-stream  velocity  ue, 

•  the  displacement  thickness 


•  the  momentum  thickness 


The  temporal  development  of  8,  8*  and  6  is  shown  in  Figure  12.  For  the  steady  inflow, 
the  thicknesses  approach  to  within  2%  of  steady  sta*e  values  by  t  —  0.23  —  0.26ms.  If  the 
approach  to  steady  state  begins  with  the  passage  of  the  upstream-facing  shock,  these  times 
give  G  ~  7  which  is  roughly  double  the  value  expected  for  a  laminar  bounda^v  layer  in  shock- 
tube  flow  [21].  This  result  seems  to  indicate  that  the  incident  shock  structure,  including  the 


contact  surface  and  the  upstream  facing  shock,  interferes  with  the  establishment  processes 
and  delays  the  approach  to  steady  state. 

The  steady  value  of  6  =  1.49mm  is  close  to  the  value  of  1.45mm  obtained  with  the 
“reference  temperature”  relations  of  Eckert  [25]  and  Equation  (7.53)  in  White  [26] .  This 
agreement  provides  some  confidence  that  the  viscous  effects  have  been  adequately  resolved. 
Also,  the  flat  plate  transition  data  from  He  &  Morgan  [27]  indicate  that  transition  to  a 
turbulent  boundary  layer  would  be  expected  after  a  distance  of  approximately  0.25m. 

For  the  transient  inflow  cases,  the  approach  to  steady  state  is  qualitatively  different  and, 
at  late  times,  there  is  a  slight  growth  in  each  of  the  thicknesses  as  a  consequence  of  the 
relaxation  of  the  free-stream  pressure  (and  hence  density). 

4.3.  Recirculation  Region 

The  form  of  the  velocity  field  in  the  recirculation  region  at  the  base  of  the  strut  appeared 
to  be  insensitive  to  the  transient  inflow.  Figure  13  shows  the  velocity  field  in  the  recirculation 
region  for  Cases  2  and  3.  The  symmetry  plane  boundary  condition  in  Case  3  allows  a  stronger 
vortex  to  form  closer  to  the  strut  base.  The  pressures  aiong  the  y  —■  0  line  in  the  recirculation 
region  are  shown  in  Figure  14  at  selected  times. 

The  length  of  the  recirculation  region,  as  measured  by  the  position  of  the  stagnation  point 
near  the  y  —  0  boundary,  is  shown  i1  Figure  15.  At  late  times,  both  the  steady  inflow  and  the 
relaxing  inflow  traces  reach  steady  values.  The  steady  state  lengths  are  Lr  -  9.5,8.3,10.1 
and  10.8mm  for  Cases  1,2,3  and  6  respectively,  which  means  that  (a)  the  symmetry  plane 
cases  have  longer  recirculation  zones  than  the  solid  wali  cases  and  (b)  the  transient  inflow 
casec  have  slightly  shorter  lengths  than  their  steady  inflow  counterparts.  Based  on  LT  for 
steady  inflow,  the  base  flow  establishes  in  a  dimensionless  time  G  ~  146  ( t  —  0.419ms)  for 
the  sol'd  wall  (Case  1)  and  G  ~  83  ( t  —  0.283ms)  for  the  symmetry  plane  (Case  6). 

Also  shown  in  Figure  15  is  the  length  of  the  recirculation  region  along  the  line  y  = 
0.0025m  which  is  half  way  up  the  base  of  the  strut.  This  length  was  determined  by  computing 
the  vertical  mass  flux 

jhy  =  J  fjuvdx,  (17) 

and  finding  the  x-location  where  mu  became  zero.  This  measure  showed  the  same  trends 
as  the  stagnation  point  measurement  except  that  G  ~  64  (f  =  0.229ms)  for  the  symmetry 
plane  (Case  6). 

4.4.  Recompression  Shock 

Figure  16  shows  the  pressure  distributions  along  the  y  =  0  plane  for  Cases  1  and 
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The  recompression  shock  is  rather  smeared  and  appears  as  the  gradual  pressure  rise  from 
x  ~  0J2m  to  x  ~  0.13m.  Observe  that,  the  pressure  distributions  for  the  transient  inflow 
(Case  2)  collapse  onto  the  one  curve  when  they  are  normalized  using  the  procedure  discussed 
in  Section  3.5. 

The  shock  position  along  a  line  through  the  midpoint  of  the  strut  base  (y  =  0.0025m)  is 
plotted  against  time  in  Figure  17.  The  shock  position  was  taken  to  be  the  grid  location  at 
which  the  pressure  rises  to  25%  of  the  inlet  pressure.  Based  on  the  length  of  the  recirculation 
zone  the  nondimensional  establishment  time  G  ~  140  ( t  ~  0.402ms)  and  G  ~  80  (f  = 
0.273ms)  for  Cases  1  and  6  respectively. 

4.5.  Wall  Loads 


The  shear  stress  and  heat,  flux  were  monitored  at  selected  points  along  the  y  =  0.02357m 
'.vail.  The  shear  stress  coefficient  is  given  by 

I 


p  ^  I  dV  wall 

W  -  T ITT’ 

and  the  heat  flux  by  the  Stanton  number 


2  Pin^in 


8T\ 

aV  \wall 


(18) 


PitxV'in Cp ( Taw  ~TW)’ 


(19) 


where  the  normal  derivatives  of  the  velocity  and  temperature  were  approximated  by  differ¬ 
ences  at  the  first  node  off  the  wall.  The  viscosity  was  evaluated  at  the  wall  temperature  using 
Sutherland’s  formula  and  pxn  and  uxn  are  inlet  conditions  computed  with  the  time  delay  spec¬ 
ified  in  Equation  (12).  Thermal  conductivity  k  was  obtained  assuming  Cp  =  1004 J/kg/K 
and  Prandtl  number  Pr  —  pCp/k  =  0.71.  The  adiabatic  wall  temperature  was  obtained 
from 


T  =  T 

-*•  aiu  *  tn 


1  +  VPr 


(7  ~  1) 


(20) 


The  development  of  these  parameters  at  x  —  0.114m  and  x  —  0.5m.  is  shown  in  Figures  18 
and  19  respectively.  Results  for  the  steady  inflow  cases  exhibit  an  approach  to  steady  state 
at  both  locations.  The  normalization  of  the  transient  inflow  measurements  is  not  entirely 
successful  as  the  values  of  both  Cf  and  Ch  settle  to  a  constant  rate  of  growth  after  some 
start-up  period. 

Normalized  pressure  histories  at  x  =  0.114m  and  0.50m  wall  are  shown  in  Figures  18(c) 
and  19(c)  for  completeness.  There  is  only  a  slight  variation  in  the  traces  for  the  two  transient 
inflow  cases  indicating  that  the  normalization  works  well  for  pressures.  The  slight  variation 
may  be  caused  by  the  use  of  the  constant  value  Uc  =  3670ma_1  in  Equation  (12). 
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4.6.  (x,t.)-Diagiams 


The  time  for  wall  pressure,  shear  stress  and  heat  transfer  to  reach  steady  state  at  a  number 
of  x-locations  along  the  y  =  0.02357 m  wall  is  shown  in  Figures  20,  21  and  22  respectively. 
An  average  of  each  quantity  between  times  t  —  Q-7?ns  and  O.Sms  was  used  as  the  Bteady 
value  and  the  settling  time  te  determined  such  that  the  flow  quantity  remains  within  2%  of 
the  steady  value  for  ie  <  t  <  0.7ms.  The  establishment  time  is  then  given  as  re  —  te  —  tuj. 
where  tu{ ,  is  the  arrival  time  of  the  upstream-facing  shock.  Only  the  steady  inflow  cases  (1 
and  6)  are  considered,  since  a  consistent  establishment  criterion  was  not  available  for  the 
transient  inflow  cases. 

The  wall  pressure  generally  settles  more  quickly  than  shear  stress  or  heat  transfer  and 
has  a  dimensionless  establishment  time  of  G  =  TeUc/x  ~  0.9  for  x  >  0.2 m.  However,  the 
establishment  is  significantly  delayed  where  the  relatively  weak  boundary-layer-interaction 
shocks  strike  the  wall  at  x  ~  0.09m.  and  0.18m.  The  pressure  measurement  does  not  settle 
at  all  where  the  relatively  strong  recompression  shock  strikes  the  wall  at  x  ~  0.3m. 

The  establishment  times  for  the  shear  stress  are  nearly  the  same  as  those  for  the  heat 
transfer  and,  for  x  >  0.4m,  can  be  approximated  by  G  ~  3.  A  similar  delay  after  the  passage 
of  the  starting  pulse  was  observed  in  the  experimental  measurements  of  East  et  al.  [24].  For 
x  <  0.3,  the  delay  beyond  the  flat  plate  value  of  G  ~  3  is  significant. 

5.  CONCLUSION 

In  order  to  study  the  establishment  of  the  major  flow  features  in  a  generic  scramjet 
combustor,  several  numerical  simulations  were  analysed.  Of  particular  interest  was  the 
effect  of  a  transient  relaxing  inflow  condition  as  found  in  a  reflected-shock  tunnel  which  has 
been  operated  in  an  undertailored  mode. 

From  the  numerical  point  of  view,  the  simulations  were  particularly  demanding  and  much 
of  the  interesting  phenomena  (such  as  turbulent  mixing  and  chemistry)  was  omitted  from 
these  simulations.  These  issues  will  be  addressed  in  a  future  study.  The  finite-difference 
scheme  generally  performed  well  but  encountered  difficulty  with  the  very  strong  shock  inter¬ 
actions  and  tended  to  produce  oscillatory  solutions  at  early  times.  With  an  incident-shock 
Mach  number  of  approximately  10,  these  oscillations  were  not  unexpected.  An  “upwinding” 
scheme  (available  in  another  version  of  the  code)  may  better  cope  with  the  strong  shocks. 

From  the  experimental  point  of  view,  it  appears  that  the  correlations  available  in  the 
literature  are  adequate  for  predicting  the  establishment  times  for  the  flow  features  of  the 
model  if  the  inflow  conditions  are  steady  and  allowance  is  made  for  the  passage  of  the  starting 
pulse  [24],  For  the  relaxing  flow,  there  are  perturbations  outside  the  2%  error  band  for  a 
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significantly  longer  time  but  these  perturbations  are  probably  smaller  than  the  measurement 
uncertainty  in  the  shock  tunnel  experiment.  Hence,  relaxing  test  flow  with  a  duration  of  the 
order  of  a  few  hundred  microseconds  is  probably  adequate  for  establishing  flow  in  moderate 
sized  scramjet  models  if  the  recirculating  zones  are  not  too  large.  The  establishment  time 
for  turbulent  jet  mixing  in  the  presence  of  combustion  needs  to  be  examined. 

The  effect  of  the  relaxing  inflow  conditions  of  the  undertailored  refiected-shock  tunnel 
seem  not  to  affect  the  overall  pressure  fields  within  the  model  nor  the  length  of  the  recircu¬ 
lation  region  to  any  great  extent.  However,  the  relaxing  flow  does  affect  features  such  as  the 
heat  transfer  and  wall  shear  stress.  It  may  also  be  expected  to  affect  the  chemical  reaction 
rates  which  will  scale  as  pLc  and  p2Lc  for  two-body  and  three-body  reactions  respectively. 
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Table  1:  Cases  Considered  in  this  Study 


Case 

Geometry 

Inflow  for  t  >  0 

Spatial  Differencing 

y  =  0  boundary 

steady /transient 

1 

short 

cold  wall 

steady 

2nd  order 

H 

short 

cold  wall 

transient 

2nd  order 

n 

short 

symmetry  plane 

transient 

2nd  order 

■ 

long 

symmetry  plane 

steady 

4th  order 

5 

long 

symmetry  plane 

transient 

4th  order 

6 

short 

symmetry  plane 

steady 

2nd  order 

Table  2:  Grid  Generation  Parameters  for  Cases  1.  2,  3,  and  6 


(a)  y-grid 


Zone 

Va 

(m) 

■©  c 

a 

P 

N 

Al/min 

(m) 

1 

0.005 

0.0 

0.5 

HR 

34 

50  x  10  -6 

2 

0.02357 

0.005 

0.5 

51 

50  x  10-6 

(b)  cc-grid 


Zone 

Xu 

(m) 

Xb 

(m) 

a 

P 

N 

AXmjn 

(m) 

1 

0.1 

0.0 

B 

1 

40 

0.001 

2 

|p|S|l 

0.1 

i 

22 

0.00025 

3 

0.124 

equal  spacing 

40 

0.00025 

4 

0.168 

0.0 

1.042 

30 

0.00025 

5 

0.168 

equal  spacing 

111 

0.00299 
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Table  3:  Grid  Generation  Parameters  for  Cases  4  and  5 


(a)  y-grid  for  x  >  —0.10964m 


Zone 

Va 

(m) 

Vb 

(m) 

Q 

0 

N 

^2/min 

(m) 

1 

0.02357 

0.0050 

0.5 

flral 

i 

51  x  10-6 

2 

0.0050 

0.0 

0.5 

Ed 

m 

53  x  10"6 

(b)  y-grid  for  —0.138m  <  x  <  —0.10964m  and  x  =  (x  +  0.138)/0. 02836 


Zone 

Va 

Vb 

a 

0 

N 

^2/min 

(m) 

(m) 

(m) 

1 

0.02357 

0.0050  x 

0.5 

1  J 

i 

65  x  10-6 

2 

0.0050  x 

-0.00639  (1  -  x) 

0.5 

O 

67  x  10-6 

(c)  x-grid 


Zone 

EH 

wGm 

Vo 

0 

N 

(m) 

1 

mm 

BUB 

8I1B1M 

1.5 

30 

0.00173 

2 

mBWlI 

3.5 

no 

0.00093 

3 

MEM 

7.0 

80 

0.00017 

4 

0.1467 

AXi  _ 

Ax;_i  ~ 

19 

0.00277 

5 

equal  spacing 

89 

0.00332 
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Table  4:  Incident  Shock  (is)  Speed  and  Upstream-Facing  Shock  (ufs)  Speed  as  defined  by 

Equation  (14). 


Case 

Vi. 

(m/ms) 

U, 

(ms) 

Vufi 

(m/ms) 

tu/» 

(ms) 

1 

3.68 

0 

2.83 

0 

2 

2.96 

0 

- 

- 

3 

2.97 

0 

- 

- 

4 

3.72 

0.0350 

- 

- 

5 

2.97 

0.0466 

- 

- 

6 

3.71 

0 

2.83 

0 
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(b) 


* 


t 


Unsteady 

expansion 


,  Test 
i  region 


Nozzle 
starting 
->-  x  process 


Figure  1:  Free-pistou  driven  reflected-shock  tunnel. 

(a)  Facility  schematic. 

(b)  Distance-time  wave  diagram  after  primary  diaphragm  rupture. 


t,  ms 


(b)  Nozzle  supply  pressure,  Ps 

100 - 

on  - 


)  I  I  _ i _ i - - - 1 - 1 

2.0  2.5  3.0  3.5  4.0  4.5 

t,  ms 


Figure  2:  Typical  pressure  histories  for  the  T4  shock  tunnel. 

(a)  Pitot  pressure  at  the  exit  plane  of  a  Mach  5  nozzle. 

(b)  Nozzle  supply  pressure. 
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Figure  4:  Typical  static  pressure  in  tbe  HYPULSE  facility. 
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Test  flow 


0.02357m 


0.1  nv 


0.138m  I  0. 1 1 4m 


0.0050m 


Figure  5:  A  generic  scramjet  combustor  with  centrally  located  fuel  injector  strut. 


(a)  Short  injector  configuration 


Supersonic 
inflow  cz; 


Supersonic  J 
inflow  psr^;;. 


i  ♦  nn 


.02357m 

WW/Z1  I 

005m  < _ _  _ 


Supersonic 

outflow 


Cold  wall 


(b)  Long  injector  configuration 

fx 

■0.138m— 0  1 14m— H 


- 0.5m 


Figure  6:  Combustor  geometry  used  in  the  computations. 
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Figure  7:  Grids  used  for  the  calculations.  Note  that  the  y- scale  is  magnified  by  a  factor  of 
five. 
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Figure  8:  Transient  inlet  conditions  as  specified  by  equations  (7)  -  (10). 
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Figure  10:  Evolution  of  the  pressure  field  for  Case  2  (transient  inflow)  at  selected  times. 
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Figure  14:  Pressure  distributions  along  y  —  0  near  the  base  of  the  injector  for  Cases  2  and 
3. 
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Figure  15:  Length  of  the  recirculation  region. 

(a)  Lr  identified  by  the  stagnation  point  location. 

(b)  length  identified  by  a  zero-croBsing  of  the  vertical  mass  flux  along  the  line  y  =  0.0025m. 
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Figure  18:  Comparison  of  loads  on  the  y  —  0.02357: 
transient  inflow. 


Normalized  pressure 


Figure  20:  Establishment  time  for  normalized  pressure  along  the  y  =  0.02357 m  wall. 
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Figure  22:  Establishment  time  for  shear  stress  coefficient  along  the  y  =  0.02357m  wall. 
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